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Abstract 

Magnetic impurity effects on metallic carbon nanotubes are studied theoretically. The resol- 
vent method for the multi channel Kondo effect is applied to the band structure of the k ■ p 
perturbation hamiltonian in the limit of the infinite onsite repulsion at the impurity site. We 
discuss the local non-Fermi liquid behavior at temperatures lower than the Kondo temperature 
7k. The density of states of localized electron has a singularity ~ M 1 / 2 which gives rise to a 
pseudo gap at the Kondo resonance in low temperatures. The temperature dependence of the 
electronic resistivity is predicted as T 1//2 , and the imaginary part of dynamical susceptibilities 
has the |cu| 1//2 dependence. Possible experimental observations are discussed. 
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I. Introduction 



Recently, carbon nanotubes with cylindrical graphite structures have been intensively investi- 
gated. Many interesting experimental as well as theoretical researches have been performed 
(see reviews [1,2] for example), and the fundamental metallic and semiconducting behaviors of 
single wall nanotubes predicted by theories [3-8] have been clarified in tunneling spectroscopy 
experiments [9,10]. 

In the magnetoresistance study of metallic nanotubes [11], disorder gives rise to the positive 
differential resistance at low temperatures due to the weak localization effects [12]. Some 
samples are known to show the negative resistance around the weak gate voltage region, which 
is interpreted possibly due to the presence of the Kondo resonance as the result of the magnetic 
impurity effect. Therefore, a magnetic impurity might change the electronic structures of 
metallic carbon nanotubes. 

The second candidate of the magnetic impurity is the fluorine (F) adatoms on the nanotubes. 
The recent experiment [13,14] of the nanographite doped with F shows appearance of electronic 
spin during the fluorine concentration, [F]/[C] = ~ 1.2. The computer simulation of F-doping 
effects on nanographite [15] supports the presence of a local spin at the fluorine atom bonded 
to the carbon atom. 

Thirdly, it is known that pristine carbon nanotubes are found to coexist with amorphous 
carbon soots [1,2]. Their effects are not of prime interests usually. However, we are aware 
that amorphous metals are modeled well by the two level system [16,17]. The low temperature 
behavior of resistivity in the two level system is similar to that of the Kondo model of the 
magnetic impurity [18,19,20]. So, the metallic nanotube with amorphous carbon soots may 
show behaviors like effects of magnetic impurities at low temperatures. This viewpoint is very 
interesting as a possible realization of the results of low temperature properties of the two level 
system, or equivalently, magnetic impurity effects. 

In this paper, we will study effects of a magnetic impurity in metallic carbon nanotubes at 
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low temperatures. The nonmagnetic impurity effects [21,22] have been studied by the present 
author, too. However, magnetic impurity effects are attractive also, when we look at our 
experience of the research of /-electron systems [23]. There are two channels of electronic 
states at the Fermi energy in metallic carbon nanotubes. As known in the magnetic systems, 
such as, heavy fermion systems [24], the non- Fermi liquid behaviors, i.e., the singular density of 
states and the power law temperature dependence of the electric resistivity, have been observed 
experimentally and explained theoretically by using the Kondo model or the Anderson model 
with multi channel scatterings. The similar effects can occur in the carbon nanotubes when 
the presence of the two scattering channels plays an important role. 

We will use the k ■ p method [25,26] for the electronic states of carbon nanotubes. The 
method well describes electronic states around the Fermi energy. The valence and conduction 
band states have the linear dispersion at the Fermi energy in metallic nanotubes. It is assumed 
that there is one magnetic impurity at a carbon site. The magnetic impurity is modeled by the 
Anderson model where there is one localized electronic state and the strong onsite repulsion is 
assumed. For the treatment of the Kondo effect, we use the resolvent method [27,28] to discuss 
the infinite repulsion case. This is sufficient for the discussion of the low temperature and low 
energy behaviors. We will solve spectral functions of the resolvents, and derive analytic formulas 
of electronic density of states, resistivity, and dynamical susceptibilities. We will discuss the 
local non-Fermi liquid behavior which could be observed at low temperatures. The density of 
states of localized electron has a singularity ~ M 1 / 2 - This singular behavior gives rise to a 
pseudo gap at the Kondo resonance in low temperatures. The temperature dependence of the 
electronic resistivity is predicted as T 1//2 , and the imaginary part of dynamical susceptibilities 
has the | cc^ | 1 / 2 dependence. Possible experimental observations are discussed. 

This paper is organized as follows. In the next section, we introduce the model and explain 
theoretical formulations. In Sec. Ill, we report the low temperature solution of resolvents. In 
Sec. IV, electronic density of states is derived. In Sees. V and VI, resistivity and dynamical 
susceptibilities are explained. Section VII is devoted to discussion, and the paper is closed with 
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summary in Sec. VIII. 



II. Model 



We will study the metallic carbon nanotubes with one Anderson impurity at the A or B sub- 
lattice site. In the total hamiltonian, 



(1) 



H tuhe is the electronic states of the carbon nanotubes, and the model based on the k • p 
approximation [25,26] represents electronic systems on the continuum medium. The second 
term ifi mp is the Anderson impurity where the onsite Coulomb interaction strength is taken as 
infinite in the course of calculation, and the resolvent formalism [27,28] is introduced in order 
to take into account of the infinite repulsion. 

The hamiltonian by the kp approximation [25,26] in the secondly quantized representation 
has the following form: 



#tube = £*fe,^fc*fc, CT > 

k,cr 

where E^ is an energy matrix: 
/ 



(2) 



E 



k 














V 



~/(k x + iky) 




(3) 



^{k x - iky) 

^(k x + iky) 


^(k x - ik y ) 

k = {k Xl k y ), and ^ a is an annihilation operator with four components. In the operator \1>, the 
fist and second columns indicate an electron at the A and B sublattice points around the Fermi 
point K of the graphite, respectively. The third and fourth columns are an electron at the A 
and B sublattices around the Fermi point K'. The quantity 7 is defined as 7 = (\/3/2)a7o, 
where a is the bond length of the graphite plane and 70 (— 2.7 eV) is the resonance integral 
between neighboring carbon atoms. When the above matrix is diagonalized, we obtain the 



dispersion relation E± = ±7^^ + where /c x is parallel with the axis of the nanotube, 

Kv4>(n) = (27r / L)(n + 0— ^/3), L is the circumference length of the nanotube, n (= 0, ±1, ±2, 
...) is the index of bands, <fi is the magnetic flux in units of the flux quantum, and v (= 0, 
1, or 2) specifies the boundary condition in the y-direction. The metallic and semiconducting 
nanotubes are characterized by v = and v — 1 (or 2), respectively Hereafter, we consider 
the case = and the metallic nanotubes v — 0. 
The second term in Eq. (1) is the impurity model: 

#im P = E d J2 dld a + y E n d,*n d - a 



+ V [X y\dUe i ^ k+K > r ^i ) +e iv e i(k + K'yr y(3) ) + hx . ]? (4) 

K,a 

where d is an operator of the localized electron with spin a; its site energy is E d ; U d is the 
onsite repulsive interaction and n d a = d\d a \ V is the mixing interaction between the electrons 
on the nanotube and the localized orbital; N s is the total number of lattice sites; r] is the angle 
between the chiral vector of the present nanotube and the tube axis direction of the armchair- 
type nanotubes; and (i — 1 — 4) is the ith component of the operator ^ cr , assuming that 
the Anderson impurity is located at the site r of the A sublattice set. When the impurity is 
located at one of the B sublattice sites, the second and fourth components of a appear in 
ifimp. However, the result of the single impurity system is the same altogether. 
The propagator of the electrons on the nanotube is defined in the matrix form: 

G(k,T) = -(T T * kta (T)^ ka (0)), (5) 

where T T is the time ordering operator with respect to the imaginary time r and ^^^(r) = 
exp(Hr)^ CT exp(— Ht). The Fourier transform of G is calculated as: 

( G# \ 



G (k,iuj n ) 



\ g k> J 



(6) 
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where u n = (2n + l)nT is the odd Matsubara frequency for fermions. The components of G 
are written explicitly: 

/ • - \ 



G K 1 (k,iuj n ) 



lUJr. 



--f(k x - ik y ) 
\ ~l(k x + ik y ) iu n 



(7) 



and 



G K )(k,iu> n ) = 



ILOr, 



-l{k x + iky) 



J 



K ~l{k x - ik y ) iu n j 



(8) 



There are several theoretical characteristic parameters. We will explain them as follows: 



• The total carbon number N s is given by 



a/3 

iV, = A x L t (—a 2 ) x 2 



4AL 



(9) 



where A is the length of the nanotube, and (v / 3/2)a 2 is the area of the unit cell. There 
are two carbons in one unit cell, so the factor 2 is multiplied. 



• The density of states near the Fermi energy E = is calculated as 



p(E) = £- [°° dk x S(E-jk x ) 



aN s 
47rZ/7o 



(10) 



Because two sites in the discrete model correspond to one site in the continuum k ■ p 
model, the density of sites in the continuum model is given by: 



p(E) 



(11) 



In the theory of Kondo effect, it is necessary to define the band cutoff D in order that 
p(E = 0) = 1/2D. Thus, the quantity D is written explicitly 



D = — 7o . 

a 



(12) 



Table I shows several combinations of parameters, D and E 1 = 7/tQ (n = 1) = (\^3ira/ L)to, 
as a reference for realistic parameter values of nanotubes. We find that D is a fairly large 



cutoff for typical single wall nanotubes. The energy E\ at the bottom of the conduction 
band n — 1 is apparently larger than the energy scale which corresponds to room tem- 
peratures or low temperatures. This would mean that effects of the upper or lower bands 
which do not cross the Fermi energy are small enough to neglect their contributions to 
the multi channel Kondo behavior in low temperatures. 

In the resolvent method [27,28], we mainly consider the limit U d — > oo so that double 
occupancy at the localized level is forbidden. The empty state and the singly occupied state 
at the magnetic impurity are taken into account in the formalism. The spin degeneracy N 
and the number of scattering channels M of the band states are generalized to take an arbitral 
integer number. In the limit of M, N, — > oo with 7 = M/N fixed, the single impurity problem 
is represented by the coupled integral equation: 

£ d M = ^- I def{e)<S> b {e + u), (13) 

7T J 

U b (u) = - fdef(e)<5> d (e + uj), (14) 

7T J 

where f = irpNV 2 , /(e) = l/[exp(e/T) + 1], 

M^n) = ~ 1 T . s (15) 

lU) n - E d - h d (lL} n ) 

is the resolvent for the singly occupied d-state, and 

*b(i"n) = ~. ^7-^ (16) 

is the resolvent for the empty rf-state. Here, we denote the even Matsubara frequency as v n . 

The validity of the multi channel problem for the application of metallic nanotube can be 
understood easily by looking at the formula of the self energies. In the case U d = 0, the self 
energy of the ci-electron is 

Z d (iu; n ) = V 2 ^-J2[ G K(k,tuj n ) + G^(k,tu n )}, (17) 
iVs k 



-8- 



where G$ is the first diagonal component of the propagator Eq. (6) and G$, is the third 
component. After inserting their explicit forms, we obtain 



= v 2 —y{ \ + 



- M ^ 2 ^-E-^ — » ( 18 ) 

where the scattering channel number is expressed as M and — D < < D with the density of 
states p = 1/2D (D is defined by Eq. (12)). Only the two bands which cross the Fermi energy 
are retained in the sum of the last line, because the characteristic energy E\ is fairly large as 
discussed above. Similarly, in the limit U d — > oo, the self energy of the resolvent <3> d is 

N * k ^' %UJn ' ~ €k 

The difference from the Ud — case is that the resolvent <!>{, appears in the right hand side. 
After taking the frequency sum, we obtain the integral form Eq. (13). Therefore, we have 
shown that the set of integral equations of the multi channel Kondo problem is valid for an 
magnetic impurity in metallic carbon nanotubes. 



III. Low temperature solution of resolvents 

The set of integral equation Eqs. (13,14) can be solved analytically in the limit of low frequency 
and low temperatures [27,28]. Before discussing physical quantities, we derive analytic forms 
of resolvents. As the resolvents mean singly occupied states and empty state physically, the 
analytic property is not the same as that of the usual propagators. We define spectral functions 
for the empty states at uj > Eq, 

A^{u)6(u - E ) = --lm^ dtb (uj + iS), (20) 

71 

and spectral functions for the occupied states at uj < E , 

Aft ee lim^oh-Im^u; + t5)}exp[f3(E - uj)], (21) 

71 
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where E is the ground state energy of the magnetic impurity at zero temperature, 5 is positive 
infinitesimal, and (5 = 1/T. 

After some calculations following Ref. [27], we find the following formula at low frequency: 

^Menr (22) 

and 

^ ~ le^r 1 . (23) 

Here, 

8M - K^)(^)]* (24) 
with 7 = M/N, and T K is the Kondo temperature: 

T K = D(^ye X p£S*). (25) 

For the special case of interests of metallic nanotubes, M = 2 and N = 2. Therefore, we 
find the singular frequency dependence around the ground state energy E : 

Af\u) ~ A^\u) ~ |£b - (26) 

This singular dependence is the main conclusion of this section. The frequency dependence is 
shown schematically in Fig. 1. 

IV. Density of states 

The local density of states of ci-electron is calculated by the convolution of spectral functions 
Ad\b ( u ) [27,28] . The explicit form of the density of states for spin a and the channel a at T = 
becomes 

p^(u>,0) * [j^^Nf][l + 6(u)U(u)+e(-u)f4u)], (27) 

where 

f±(u) = a±\u\ A °» + b±\u\ A °\ (28) 



with 




7r)sin(7rA ch ) J B(2A sp ,A ch ), 



(29) 



— cos(7rA ch )a_, 



(30) 



b. 



+ 




7r)sin(7rA ch )5(2A ch ,A sp ), 



(31) 



cos(7rA ch )6 + . 



(32) 



Here, u 



[(1 + ^)/^](uj/T k ), B(x,y) is the Beta function, and = irT K /F is the weight of 



channel fluctuations. Further, A sp = 1/(1+7) an d A ch = 7/(1+7) are the scaling dimensions of 
spin and channel fields, respectively. As A ch oc M, this measures the magnitude of fluctuations 
from the channel degree of freedom. Also, A sp oc N means that this is a measure of contributions 
from the spin degree of freedom. Both quantities determine the degree of singularities of 
electronic density of states and physical quantities at low frequencies. They are the most 
important parameters introduced in this section. 

Specially for metallic carbon nanotubes, we know that A sp = A ch = 1/2. This implies the 
singularity around the Fermi energy uj = 0: 



Such the singular functional form implies that a pseudo gap develops at the top of the Kondo 
resonance peak which appear at temperatures much lower than T K . There appears a dip in 
the density of states at the Fermi energy. The dip structure of the density of states is shown 
schematically in Fig. 2. This is the local non-Fermi liquid behavior discussed in detail in the 
literatures [24,29]. If it is possible to measure the local density of states of a metallic atom 
attached to the carbon nanotubes, for example, by scanning tunneling microscope, we could 
observe such the pseudo gap behavior when role of the multi channel scatterings is dominant. 



1 + 6(u)\u\* + 0(-u)\u\*, 




(33) 
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V. Resistivity 

In this section, we consider the electric resistivity in order to look at how the singular behavior 
will be observed. The scattering rate r is calculated from the scattering t matrix: 



(34) 



r^T)- 1 = -2ImtW(u + i8,T), 

2Tp^ a (uj,T) 
pN 

The relation with the electronic resistivity 

m ~ [J de(-^)T(e,T)}-i (35) 
gives the low temperature behavior: 

© ~ 1 - c(^r in(Asp ' Ach) + (36) 
p(0) T K 

where c is a constant, but it is difficult to obtain its explicit form only from the information of 
^-dependence of p u>a . Here, we stress that the above singular functional form agrees with that 
obtained from the conformal field approach [30,31]. Therefore, the temperature dependence of 
the leading term (T/Tx) mm< - Asp ' Ach - ) is a general result which is independent from method of 
theoretical treatment. 

For the metallic carbon nanotubes, we already know A sp = A ch = 1/2. Therefore, the low 
temperature behavior 



is expected from the above general formula. 

VI. Magnetic susceptibility 

Spin and channel susceptibilities are calculated by the linear response function. The spin 
susceptibility is the magnetic susceptibility in other words. The imaginary parts of dynamical 
susceptibilities are defined as 

Xs P = Jj lm Xs P (38) 



and 



Xch = ^ Im Xch- (39) 
The first term of x'L at T = is calculated to be 



a 



X;(^,0)~^sgnu;p|( A -- A -), (40) 

J- k 



where 

C SP = 7 A s 2 p sin(7rA sp ) J B(A sp , A sp ). (41) 
The second correction gives the u dependence: 

X>; p (uj,0)~\^-^\ (42) 
In the similar way, the dominant term of x"h at T = becomes 

^0)~^sgncc#|( A -- A -\ (43) 

J-K 

where 

C ch = W^A sp sin(7rA ch )5(A ch , A ch ). (44) 

The second term has the u dependence: 

&(u>,0) ~ \u,\V**-*"\ (45) 

The above general formulas Eqs. (40-45) reduce to that of metallic carbon nanotubes. The 
result of singular behavior is common for spin and channel susceptibilities: 



x"(uj, 0) ~ Asgnuj(l - Bsj M + ...), (46) 

where A and B are constants. We find \J\u\ dependence at low frequencies. We note that the 
term sgn oj is related with the real part: 
-// m\ i r max(cJ,T) 1 

X \u, T) ~ -log[ ^-^]. (47) 



VII. Discussion 



We have discussed the solutions of the integral equations in the limit M, N — > oo in this paper. 
On the other hand, the higher order corrections of 0(1/ N) [27,28] are known to change the 
coefficients of the leading terms of singular behaviors, but they do not affect the power of 
singularities. Therefore, the non-Fermi liquid behaviors are not affected by the higher order 
corrections by the 0(1 /N) expansion. 

In the course of the investigation, we have neglected the valence and conduction band 
states which do not cross the Fermi energy. This approximation applies well at low energies 
\uj\ ~ T K <C Ei. Typical magnitudes of E 1 have been shown in Table I. However, at higher 
energies or at high temperatures, such the deep bands will give some corrections to low energy 
singular behaviors. Furthermore, functional forms of spectral functions at entire frequencies 
can be affected by the presence of deep bands. In order to treat the deep band effects, we have 
to solve the resolvent equations numerically at least, even though it is beyond the scope of this 
paper. 

There are two scattering channels when the impurity interacts with the single wall metallic 
nanotube. Is it possible that there are more scattering channels? When two metallic carbon 
nanotubes are present and one impurity interacts with both nanotubes, four scattering channels 
are realized as long as the interactions between electrons of two nanotubes are negligible. The 
powers of the singularities of the density of states and physical quantities become different from 
those of the case of the single nanotube. By putting M = 4 and N = 2, the general formulas of 
this paper can be applied to one magnetic impurity which interacts with two aligned metallic 
nanotubes. If the impurity interacts with more nanotubes, such the general formulas are useful 
in order to predict singular behaviors which depend on the scattering channel number. 

In the literature, for example in [32], the power law temperature dependences of electronic re- 
sistivity have been interpreted by the Luttinger liquid picture of the correlated one-dimensional 
systems. The Kondo effect in the Luttinger liquid has been studied by Furusaki and Nagaosa 
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[33]. A further study of the Kondo impurity with two scattering channels in the presence of 
conduction bands with Luttinger liquid properties is a fascinating extention of the present work. 

VIII. Summary 

Magnetic impurity effects on metallic carbon nanotubes have been investigated theoretically. 
The resolvent method for the multi channel Kondo effect has been applied to the band structure 
of the k ■ p model in the limit of the infinite onsite repulsion at the impurity site. We have 
discussed the local non-Fermi liquid behavior which could be observed at temperatures lower 
than the Kondo temperature Tk- The density of states of localized electron has a singularity 
~ | a; 1 1 I 2 . This singular behavior gives rise to a pseudo gap at the Kondo resonance in low 
temperatures. The temperature dependence of the electronic resistivity is predicted as T 1//2 , 
and the imaginary part of dynamical susceptibilities has the M 1 / 2 dependence. 
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TABLE I. Parameters D and Ei for typical metallic carbon nanotubes. 



Nanotube 


D (70 : 


= 2.7 eV) 




Ei 




(5,5) 


8.60 70 


= 23.2 eV 


1.99 


7o = 


5.37 eV 


(10,10) 


12.2 7o 


= 32.9 eV 


1.40 


7o = 


3.78 eV 


(15,15) 


14.9 70 


= 40.2 eV 


1.15 


7o = 


3.11 eV 


(6,0) 


16.3 70 


= 44.0 eV 


1.05 


7o = 


2.84 eV 


(9,0) 


24.5 7 o 


= 66.2 eV 


0.698 


7o = 


1.88 eV 
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Figure Captions 

Fig. 1. The singular behavior of the spectral functions A^(ou). The figure is drawn schemati- 
cally. 

Fig. 2. The dip structure of the density of states p(u). The figure is drawn schematically. 
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